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Abstract: Cellular senescence plays an important role in the age-dependent functional decline of organs and organ 
systems, as well as in age-related pathologies, such as cancer. Therefore, a better understanding of its underlying 
molecular mechanisms is crucial in the search for intervening measures. In this study, we considered the role of DNA 
methylation in senescence. We found that senescence is associated with global DNA hypomethylation, but also involves 
site-specific DNA hypo- and hypermethylation. In some cases, this differential methylation may affect gene expression and 
thereby modulate functional processes within cells. However, the majority of the CpG sites that were differentially 
methylated did not correspond with altered gene expression, suggesting that DNA methylation affects senescence by 
other means also, such as, for instance, genome stability. 



INTRODUCTION 

Cellular senescence is associated with the terminal cell 
cycle arrest of cells, and was first described by 
Hayflick and Moorhead, who observed that human 
diploid fibroblasts exhibit a limited in vitro replicative 
lifespan of approximately 50 ± 10 population 
doublings, after which the cell cultures progressively 
decline [1, 2]. Later, such senescent cells were shown 
to accumulate in tissues in vivo with increasing age [3, 
4], and when primary cells from species with different 
maximum lifespans were taken into a cell culture, their 
in vitro lifespan correlated with the maximum lifespan 
of the donor species [5, 6]. These observations 
supported a role for cellular senescence in organismal 
aging. 

There are several different ways, in which cellular 
senescence may contribute to the age-related functional 
decline of organs and organ systems and age-related 
pathologies. On the one hand, the senescence of cells 
that rely on proliferation for their proper function, such 
as stem cells and immune cells, limits the repair 
capacity in aging tissue [7] and contributes to the age- 



dependent deterioration of the immune system [8]. On 
the other hand, while the pro-inflammatory micro- 
environment of senescent cells can contribute to the 
development and progression of inflammatory diseases 
and cancer [9, 10], the failure of senescence 
mechanisms can also result in malignant transformation 
and cancer progression [11]. Further, a systems biology 
approach showed the similarity between genes involved 
in cellular senescence, aging and age-related 
pathologies, suggesting these may have co-evolved 
[12]. Thus, understanding the mechanisms that underlie 
senescence constitutes a major step towards the goal of 
ensuring healthy aging. 

While certain molecular changes have already been 
associated with the induction of senescence, such as 
telomere shortening [13, 14], and the increased 
expression of senescence-associated genes like 
CCND1 [15, 16], the role of epigenetic phenomena 
[17] and changes in chromatin structure have recently 
also been considered for their roles in aging and 
senescence (see also [18-20]). Narita et al. [21] were 
the first to describe the focal formation of 
heterochromatin in senescent cells, which they termed 
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"senescence-associated heterochromatin foci" (SAHF). 
They showed that this heterochromatin formation was 
targeted to E2F target genes by RB/E2F and reinforced 
terminal cell cycle arrest through the permanent 
downregulation of cell cycle regulators. Swanson et al. 
[22], meanwhile, described an early step in the 
establishment of senescence that they tenned the 
"senescence-associated distension of satellites" 
(SADS), which was caused by the relaxation of the 
chromatin structure in pericentromeric satellite regions. 
Thus, both the loss of heterochromatin in regions of 
constitutive heterochromatin and the gain of 
heterochromatin in previously transcriptionally active 
regions may contribute to the establishment of 
senescence. 

Several recent studies involved massive parallel 
sequencing projects that compared the DNA 
methylomes of blood cells from donors of varying 
ages ranging from newborns to centenarians [23-25]. 
They commonly observed a decrease in global DNA 
methylation with increasing age. This loss of DNA 
methylation was observed in diverse regions of the 
genome and affected various functional groups of 
genes. Heyn et al. (2012) found that mainly CpG-poor 
and tissue-specific promoters were affected by DNA 
hypomethylation, while McClay et al. (2014) observed 
reduced DNA methylation levels in regions associated 
with polycomb proteins or activating histone marks. 
DNA hypermethylation, on the other hand, occurred 
more frequently in CpG-rich sequences, and seemed to 
be more site-specific than DNA hypomethylation. In 
addition, a comparison of the DNA methylomes of 
different cell types at different passages showed that 
the differential methylation of several sites served as a 
bio marker for senescence [26]. 

Even though DNA hypomethylation seemed to be less 
site-specific than DNA hypennethylation, the 
observations that the loss of DNA methylation only 
occurs in senescing cell strains, and not in immortalized 
cells [27], and that the inhibition of DNA methylation 
can induce growth arrest in immortal cells [28], suggest 
that DNA hypomethylation does play an important role 
in the establishment of senescence. 

Here, we investigate the role of differential DNA 
methylation in senescence in WI-38 human diploid lung 
fibroblasts. Senescent WI-38 cells exhibited down- 
regulation of DNMT1 and reduced global DNA 
methylation. The extensive changes to the DNA 
methylation profile are shown to affect partially distinct 
functional groups of genes according to whether there is 
hypomethylation or hypermethylation. However, 
changes in DNA methylation only corresponded with 



differential gene expression in a small fraction of the 
affected sites. 

RESULTS 

Senescence is associated with reduced global DNA 
methylation 

To date, several studies have shown age-dependent 
alterations in the DNA methylation profile of cells. 
While both hypermethylated and hypomethylated sites 
have been detected with the increasing age of the 
subjects, DNA hypomethylation was more frequent [23- 
25]. Our previous results indicated that DNA 
methyltransferase DNMT1, which perforins 
maintenance and de novo DNA methylation and is the 
predominant DNA methyltransferase in somatic cells 
[29, 30], and Methionine adenosyltransferase 2 alpha 
(MAT2A), which plays a role in the formation of S- 
adenosylmethionine (SAM), the substrate for DNA 
methylation, were downregulated in senescent cells 
(manuscript 000141R1 accepted in AGING). This was 
confirmed at the protein level (Fig. 1A, B) and 
correlated with decreased global DNA methylation, as 
shown by the increased incorporation of radioactively 
labeled dCTP into DNA from senescent cells (Fig. 1C). 

Hypomethylation of CpG sites in the promoter 
regions of genes is observed with increasing 
senescence ratio of cultures 

In order to study how this deregulation of DNMT1 and 
the global DNA hypomethylation may affect the DNA 
methylation patterns at specific sites, genomic DNA 
was extracted from cultures of different PD levels and 
assessed for the methylation status of CpG sites with 
regulatory functions in the promoter regions of genes 
using Illumina® HumanMethylation27 BeadChips, 
which interrogate 27,578 CpG sites spanning 14,495 
genes (for detailed results, see SI). 

When comparing the PD 54 to PD 38 cultures, the array 
data showed that 1,849 CpG sites were affected by 
differential methylation (Fig. 2A); these sites were 
predominantly hypomethylated (66.7% of all 
differentially methylated sites) (Fig. 2B). This 
increasing hypomethylation of CpG sites is in line with 
the observed reduction in DNMT1 expression and 
global DNA hypomethylation in more senescent 
cultures (Fig. 1). 

When comparing the PD 47 and PD 38 cultures, the 231 
CpG sites affected by differential methylation appear to 
be randomly distributed across the genome, whereas 
when comparing PD 54 and PD 38, several clusters with 



www.impactaging.com 



565 



AGING, July 2014, Vol. 6 No.7 



£ E 



1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 



\\\ 



< 

fN 
I- 
< 



1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 



PD38 



PD47 



PD54 



PD38 



PD 47 



PD 54 



DNMT1 

MAT2A 
GAPDH 



PD38 



r- - 



PD47 



PD 54 



1.6 
' 1.4 
. 1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 



PD 38 



PD47 



PD 54 



Figure 1. Reduced protein levels of DNMT1 correlate with decreased global DNA methylation. (A) 

DNMT1 and MAT2A protein levels normalized to protein levels in PD 38 cells. Averages of three biological 
replicates ± standard deviation. (B) Western blot images. (C) Cytosine extension assay showing average 
incorporation of radioactively labeled dCTP from three samples per PD level; error bars indicate standard. 



increased incidence of hypomethylated or hyper- 
methylated CpG sites were found (Fig. 2C). Examples 
of such clusters are shown in Figure 2C (i, ii, iii). The 
observed DNA methylation changes were detected 
predominantly in regions of moderate or high gene 
density, probably due to the enrichment of probe sets 
specific to those regions on the BeadChip that was used 
for this study. Therefore, the absence of DNA 
methylation changes detected in gene-poor regions of 
the genome does not allow for any conclusions about 
changes in the DNA methylation status in those regions. 
Interestingly, in the PD 54 and PD 38 comparison, only 
42 genes that showed differential promoter methylation 
were also affected by changes in gene expression (Table 
1), while none of the genes with differentially 
methylated promoter regions were found to be 
differentially expressed when comparing the PD 47 and 
PD 38 cultures. Among those 42 genes, for 24 genes the 
change in gene expression corresponded with the 
change in DNA methylation. These genes included 
CDK2 and E2F2, both regulators of the Gl/S 
progression, which were transcriptionally repressed and 
showed hypermethylation in their promoter regions 
(Table 1). Hypermethylation of the promoters of these 



genes may mediate the permanent cell cycle arrest of 
senescent cells in the G1/G0 phase of the cell cycle. 

Further, several genes that function as transcription 
factors or regulate the activity or target specificity of 
transcription factors were among the 24 genes with 
corresponding changes in gene expression, E2F2, 
FHL1, which enhances the transcriptional activity of 
NFATlc [31], IDS, which inhibits several transcription 
factors of the basic -helix-loop-helix family [32], 
NFE2LS, which was identified as a negative regulator 
of the antioxidant-response element containing 
promoters [33], NR4A2, which can be transcriptionally 
induced in response to pro-inflammatory signals [34] 
and may then induce the expression of additional 
cytokines through the interaction with NF-icB/p65 [35], 
RELB, which is a transcription factor of the NF- 
kB/RELB family and represses pro-inflammatory gene 
expression in fibroblasts [36], and TGIF2, which 
transcriptionally represses TGF(3 target genes [37]. This 
supports the involvement of multiple transcriptional 
programs in the regulation of the senescence-associated 
gene expression profile as well as the establishment of 
the senescence-associated phenotype. 
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Table 1. Genes affected by differential expression and DNA methylation. 

Numbers indicate log2-fold changes in the transcript level or CpG methylation level when comparing PD 
54 to PD 38 cultures; information on gene function retrieved from Genecards [53]. 





Gene 


DNA 


Gene function 




expression 


methylation 




ALDH1A3 


-0.41 


0.78 


Aldehyde dehydrogenase, which binds to retinal 


APP 


0.72 


-1.11 


Amyloid beta (A4) precursor protein, best studied 
for its role in Alzheimer's disease 


ATP5D 


0.48 


1.1 


Mitochondrial membrane ATP synthase 


CDK2 


-0.43 


0.55 


Cyclin-dependent kinase, mediates Gl/S 
progression together with CCNE and the transition 
from S to G2 phase with CCNA 


COL16A1 


-0.43 


-0.68 


Type XVI Collagen, involved in mediating cell 
attachment and morphology 


CPS1 


-0.59 


-0.52 


Mitochondrial Carbamoyl Phosphate Synthase 1, 
involved in the urea cycle 


DPP4 


0.46 


0.4 


Dipeptidyl-Peptidase 4, membrane glycoprotein, 
serine exopeptidase 


E2F2 


-0.6 


1.94 


Transcription factor involved in the regulation of 
Gl/S transition and DNA replication 


FHL1 


0.57 


-0.42 


Contains zinc finger domains and LIM domains, 
may be involved in muscle development 


GFPT2 


-0.53 


2.48 


Controls the flux of glucose to the hexosamine 
pathway 


GPI 


0.64 


0.73 


Glucose-6-Phosphate Isomerase, involved in 
Glycolysis 


GST02 


-0.45 


0.49 


Glutathione S-Transferase Omega 2, may be 
involved in ascorbic acid recycling 


HMGB2 


-1.14 


-1.94 


Preferentially binds to ssDNA and may bend it, 
co-factor to RAG during V(D)J recombination 


HYLS1 


-0.7 


-0.97 


Associated with centriole, involved in cilia 
formation 


ID3 


1.02 


-1.36 


Inhibits DNA binding of transcription factors, 
including E2A 


KIF4A 


-0.72 


-0.95 


Motor protein, involved in spindle organization 


MLLT11 


0.5 


-0.99 


Fused with a number of translocation partners in 
leukemia 


MOV10 


-0.43 


-1.53 


RNA helicase, involved in RISC-mediated post- 
transcriptional gene silencing 


MVP 


0.51 


0.51 


Major vault protein, plays a role in the formation 
of scaffolds for signal transduction 


NEFH 


-0.51 


0.56 


Neurofilament protein 


NFE2L3 


-0.76 


1.5 


Transcription factor, which binds to antioxidant 
response element in target promoters 


NR4A2 


0.83 


-1.65 


Nuclear Receptor Subfamily 4 (NURR1), may 
function as a transcription factor 


NTF3 


-0.61 


0.54 


Neurotrophin 3, involved in neural survival 


PITX1 


1.43 


0.61 


Transcriptional regulator, involved in the 
regulation of development 


PPIL5 


-0.5 


1.03 


Negative regulator of TNFRSF9 signalling 


RAB11FIP5 


0.4 


-2.1 


Involved in protein trafficking 


RAM PI 


-0.61 


1.16 


Membrane protein, co-receptor 


RELB 


-0.44 


0.5 


Dimerizes with NF-kB to modify its preference for 
transcription target binding 


RELN 


-0.65 


-1.07 


Extracellular matrix serine protease 


RPL39L 


-0.62 


0.86 


Ribosomal protein L39 like 
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SCAMP3 


0.45 


1.84 


Involved in post-Golgi recycling pathways and 








protein trafficking 


SLC4A4 


0.44 


-0.59 


Sodium/bicarbonate co-transporter 


SLIT2 


-1.33 


-1.95 


May regulate cell migration 


SPATA18 


0.55 


1.23 


Involved in the repair or degradation of damaged 








mitochondria 


SSRP1 


-0.61 


0.45 


Part of the FACT complex, which is involved in 








destabilizing and reassembling nucleosomes 








during transcription, replication and DNA repair 


TFRC 


-0.51 


1.19 


Receptor involved in iron uptake 


TGIF2 


-0.52 


0.56 


Transcriptional repressor of TGFbeta responsive 








genes, through the recruitment of histone 








deacetylases 


TMEM35 


0.75 


-0.67 


Transmembrane protein 35 


TMEM47 


0.64 


-1.02 


Localized to ER 


TMSL3 


-0.45 


-0.41 


Unknown function 


TRO 


-0.49 


-0.71 


Involved in cell adhesion 


ZAK 


0.54 


1.57 


Stress-activated, involved in JNK and p38 








pathways 



While 52.3% of all differentially methylated CpG sites 
were located in CpG islands, among the 42 CpG sites 
that resided in promoters of differentially expressed 
genes a higher fraction of sites was located in CpG 
islands - 87.5% of sites which corresponded with the 
change in gene expression and 72.2% of sites which did 
not correspond with the change in gene expression. 
Additionally, differentially methylated sites that 
corresponded with a gene expression change were 
located slightly closer to the transcription start site 
(average distance of 412 nucleotides) when compared to 
all differentially methylated CpG sites (average distance 
of 421 nucleotides), while differentially methylated 
sites that did not correspond with the observed gene 
expression change were located slightly more distant 
(average distance of 495 nucleotides). However, due to 
large variations in the distance to the transcription start 
site, these differences were not statistically significant. 
Further, fifteen out of the 24 cases, in which differential 
DNA methylation corresponded with differential gene 
expression, consisted of hypermethylated CpG sites and 
corresponding repression of the transcript. 

Thus, while CpG sites located in CpG islands were 
more frequently associated with differential gene 
expression, the majority of the observed changes in 
DNA methylation did not correlate with changes in 
gene expression. These changes may not affect CpGs 
that directly regulate the gene expression, but they 
may modify the accessibility of the promoter region 
to transcriptional regulation or affect genome 
stability. 



In order to better understand the potential functional 
implications of the differential methylation of all the 
sites affected, genes with differential methylation status 
of CpG sites within their promoter regions were 
functionally classified using the DAVID and FunNet 
software. This showed that both hypo- and hyper- 
methylated CpG sites are found in the promoters of 
genes involved in cell cycle regulation, apoptosis and 
senescence, as well as in epigenetic and transcriptional 
regulation (Fig. 3A). When comparing the PD 54 and 
PD 38 cultures, the remaining genes were enriched for 
the biological processes of structural development, cell 
adhesion and localization, transmembrane transport, 
response to extracellular stimuli and cell signalling and 
maintenance of cellular homeostasis and ion 
homeostasis as determined using the g:Profiler 
software. 

When comparing the sites that were differentially 
methylated between PD 47 and PD 38 to those 
differentially methylated between PD 54 and PD 38, 
52.6% of the sites affected in the PD 47 and PD 38 
comparison were found to be common. This may 
indicate that those sites are subject to differential 
methylation early during differentiation or the 
establishment of senescence. To better understand what 
processes may be affected by differential methylation in 
both PD 47 and PD 54 cultures, detailed functions of 
the affected genes were gathered from Genecards and 
functionally classified (Fig. 3B). Among the genes with 
hypomethylated CpG sites in their promoter, 
involvement in transcriptional or epigenetic regulation, 
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Figure 2. Changes in DNA methylation pattern in 
senescent cells. (A) Number of CpG sites affected by 
differential methylation comparing PD 47 and PD 54 to PD 38 
cultures. (B) Number of CpG sites in the given comparisons 
affected by hypo- (dark grey) or hypermethylation (light grey). 
(C) Distribution of CpG methylation changes across the 
chromosomes. Top: changes between PD 47 and PD 38; Bottom: 
changes between PD 54 and PD 38. Black/grey lines indicate 
gene densities, green lines represent CpG sites hypomethylated 
in the older cultures and red lines indicate CpG sites 
hypermethylated in older cultures. 



cell signalling, cell adhesion and transport of ions were 
prevalent, whereas among the genes with hyper- 
methylated promoters, the major functional groups 
affected were cell signalling, transcriptional regulation 
and cell metabolism. This resembles the functional 
groups affected by early gene expression changes and 
supports the conclusion that transcriptional and 
epigenetic regulation may affect the establishment of 



senescence at the level of gene expression (manuscript 
000141R1 in AGING). Changes in cellular signalling 
may also be involved in the establishment of senescence 
by modifying the activation of transcriptional programs 
in response to extracellular and intracellular stimuli, 
such as oxidative stress and a pro-inflammatory 
micro environment, as well as by inducing changes in 
the cell physiology. 

DISCUSSION 

As reported (manuscript 0001 4 1R1 accepted in 
AGING), our previous study indicated that in addition to 
the numerous transcription factors, 2.7% of all the 
differentially expressed genes found when comparing 
PD 54 and PD 38 were associated with functions in 
epigenetic regulation. Among those, DNMT1 was 
downregulated in senescent cells, which correlated with 
the global reduction of DNA methylation. This supports 
the conclusion that there is a senescence-dependent loss 
of heterochromatin, as also observed in different tissues 
of aging model organisms [25, 27, 38-40]. A study by 
Fairweather et al. (1987) showed the treatment of 
human foetal lung fibroblasts (MRC5) with 5- 
Azacytidine to result in a reduction of the content of 
methylated CpGs in the genomes of the cells and 
correlate with a reduced in vitro lifespan, thus 
suggesting that DNA hypomethylation may play an 
important role in the control of the lifespan of human 
diploid fibroblasts [41]. This effect may be mediated by 
loss of DNA methylation in promoter regions resulting 
in the deregulation of gene expression profiles, as 
proposed in the heterochromatin-loss model of aging 
[42]. On the other hand, the loss of DNA methylation in 
intergenic regions and at repetitive sequences may 
result in genomic instability through a loss of silencing 
of those regions [43]. 

When we considered site-specific differences in DNA 
methylation when comparing senescent and 
proliferatively active cell cultures, we found extensive 
differential methylation of CpGs in the promoter 
regions of genes (Fig. 3). Two thirds of those sites were 
affected by hypomethylation, while one third of the 
sites were hypermethylated. These observations are in 
line with the finding that aging is associated with a 
global loss of DNA methylation [25]. Heyn et al. (2012) 
further showed that differentially methylated CpG sites 
mostly occurred in the intergenic and intronic regions, 
with only 10% affecting promoter regions [25]. As the 
Illumina J HumanMethylation27 BeadChip contains 
probes that detect methylated CpGs in promoter 
regions, our dataset may only represent a fraction of the 
differential CpG methylation associated with 
senescence in WI-38 cells. 



www.impactaging.com 



569 



AGING, July 2014, Vol. 6 No.7 



A 



Hypomethylated CpG sites 



Hypermethylated CpG sites 



PD 38- PD47 



PD38-PD54 




B Hypomethylated CpG sites Hypermethylated CpG sites 




lism 



Figure 3. Functional classification of the genes affected by differential DNA methylation. 

(A) Pie charts showing functional classification of hyper- and hypomethylated CpG sites. (B) Pie charts 
showing functional classification of CpG sites affected by methylation changes when comparing PD 47 
and P38 as well as PD 54 and PD 38 cultures. 



Since the probes on the BeadChip mainly covered CpGs 
in promoter regions, it is not surprising that only limited 
events of differential methylation were detected in 
gene-poor regions (Fig. 3). However, when comparing 
the DNA methylation patterns in PD 54 to PD 38 cells, 
some chromosomal regions with clustered events of 
CpG hypo- or hypermethylation became apparent. 

Surprisingly the overlap between the DNA methylation 
profiling and gene expression profiling datasets is very 
limited. Only 24 out of the 1,849 differentially 
methylated CpG sites corresponded with the differential 
expression of the according gene. This indicates that in 



this model, DNA methylation changes are not the major 
regulatory factor of the senescence-associated 
deregulation of transcription, which is in contrast to the 
proposed heterochromatin-loss model of aging [42]. 
However, there were several transcription factors 
among the genes for which differential DNA 
methylation correlated with differential gene 
expression, as well as CDK2 and E2F2, two central 
regulators of cell cycle progression. These genes 
exhibited promoter hypermethylation along with 
reduced transcript expression, indicating that the hard- 
wiring of the transcriptional repression through DNA 
methylation may play a role in the maintenance of the 
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permanent growth arrest. The observation that 1 5 out of 
the 24 genes were affected by promoter hyper- 
methylation and transcriptional repression, including 
CDK2, E2F2, TGIF2, RELB and NFE2L3, suggests that 
targeted promoter hypermethylation may affect the 
regulation of the gene expression profiles to a greater 
extent than the loss of DNA methylation through 
inducing a more permanent suppression of the 
expression of cell cycle and transcriptional regulators. 
This is in line with the observation that DNA 
hypermethylation seems to be more targeted than DNA 
hypomethylation [24]. 

The results of the functional classification of the genes 
affected by differential promoter methylation indicated 
that the processes affected were more diverse than those 
observed for the changes in gene expression. The 
differentially methylated promoters when comparing 
PD 54 and PD 38 showed enrichment in KEGG 
pathways involved in metabolism (arachidonic acid, 
glycine/serine/threonine, linoleic acid, nitrogen and 
glycerolipid metabolism and steroid hormone 
biosynthesis) and cell signalling (Cytokine-cytokine 
receptor interaction, NOD-like receptor signalling, 
JAK-STAT signalling, and Toll-like receptor 
signalling), while functions in the processes of cell 
cycle, DNA repair and transcriptional and epigenetic 
regulation contributed smaller fractions of differentially 
methylated CpG sites. The overrepresentation of those 
metabolic and cell signalling pathways corresponds 
with previously reported changes in lipid metabolism 
[44] and in the expression of cytokines and TLRs in 
several aging non immune system cells [45, 46]. Thus, 
promoter hypomethylation in those cases may increase 
the accessibility for transcriptional regulators without 
inducing permanent transcriptional activation, as no 
changes in transcript levels were detected. 

When comparing the sites that were affected by 
differential methylation between PD 47 and PD 38 to 
those for PD 54 and PD 38, 52.6% of the PD 47-to-PD 
38 differentially methylated CpG sites were found to be 
consistently differentially methylated. The hyper- 
methylated CpGs among the overlapping sites mainly 
affected genes with functions in transcription, 
metabolism or cell signalling, while the hypomethylated 
CpGs affected a higher variety of functional groups, 
including transcriptional regulation, cell signalling, cell 
adhesion, transporters and epigenetic regulators. 

Further, this 52.6% overlap in differential methylation 
is considerably smaller than the 80% overlap observed 
for the gene expression profiles (manuscript 000141R1 
accepted in AGING), suggesting that changes in DNA 
methylation may be more dynamic across the cellular 



lifespan. In addition, several of our findings indicate a 
role of DNA methylation in the establishment of 
senescence that does not solely depend on the regulation 
of gene expression - the majority of the differentially 
methylated CpG sites when comparing PD 54 and PD 
38 cultures were hypomethylated, while the majority of 
differentially expressed genes were downregulated, and 
only 24 out of 1,849 differentially methylated CpG sites 
corresponded with altered gene expression. Thus, it will 
be important to determine the effects of the extensive 
changes in DNA methylation profiles on genome 
integrity, nuclear architecture and nuclear matrix 
attachment, and their role in cellular senescence in 
future studies. It would also be important to analyse the 
roles of various histone modifications in senescence- 
associated gene expression regulation [47]. 

In summary, senescence in WI-38 cells was associated 
with the downregulation of DNMT1 and the global 
DNA hypomethylation. At the level of specific CpG 
sites, two-thirds of the differentially methylated sites 
were hypomethylated. However, differential DNA 
methylation did not correspond with differential gene 
expression in the majority of cases, suggesting that 
differential DNA methylation may affect the 
establishment of senescence through other mechanisms 
than transcriptional regulation. Further studies are 
needed to shed more light on the possible mechanisms 
involved and dissect the precise nature and biological 
repercussions of DNA methylation changes in 
senescence. Our datasets can be used as an important 
roadmap for the future analysis of the roles of DNA 
methylation in senescence. 

METHODS 

Cell culture. Human foetal lung fibroblasts (WI-38, 
ATCC, CCL-75TM) were maintained in HyClone 
minimum essential medium (MEM) Alpha Modification 
(ThermoScientific) containing 10% (v/v) foetal bovine 
serum (FBS) (Gibco) in a humidified Forma Steri-Cycle 
C0 2 Incubator (ThermoScientific) at 37°C and 6% C0 2 . 
In order to compare the senescence-associated changes, 
three different population doubling (PD) levels of cells 
were used: PD 38 (proliferatively active), PD 47 (pre- 
senescent) and PD 54 (senescent), as described 
previously (manuscript 000141R1 in AGING). 

Cytosine extension assay. DNA was isolated from cells 
from two 1 0-cm cell culture plates per sample and three 
samples per PD level using the DNeasy® Blood and 
Tissue Kit (Qiagen) and a cytosine extension assay was 
performed as previously described [48]. A PerkinElmer 
Liquid Scintillation Analyzer Tri-Carb 2910 TR was 
used to measure CPM values. 
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DNA methylation profiling. Genomic DNA (gDNA) 
was isolated from cells using the DNeasy® Blood and 
Tissue Kit (Qiagen). Three samples per PD level were 
used, with each sample consisting of cells harvested 
from two 10 cm cell culture plates. 500 ng of DNA was 
used to determine DNA methylation profiles on 
Illumina HumanMethylation27 BeadChips, according^ 
to the manufacturer's protocol. Illumina" 
GenomeStudio software was used for the determination 
of beta values (percentage of methylation of specific 
CpG sites). Differential DNA methylation analyses 
were performed using an Illumina Custom model, 
which produced Diff Scores as a measure of 
significance. Diff Scores of -13/13 were set as cut-off 
for significance. 

Functional classification. Functional classification of 
genes was done using several different softwares: 
FunNet Transcriptional Networks Analysis 1 , g:Profiler 
[49, 50], and DAVID Bioninformatics Resources 6.7 
[51, 52] and compared with information on the genes 
available from the Genecards database [53]. All 
differentially methylated CpG sites were counted 
individually. 

Western immunoblotting. Protein was extracted by 
sonication of cells in 100 uL cold 1% SDS containing 
protease inhibitor (Roche). 40 ug of protein was 
separated by SDS-PAGE in slab gels of 6% (DNMT1) 
or 10% (MAT2A) polyacrylamide and transferred to 
Amersham Hybond-P polyvinylidene difluoride 
(PVDF) membranes (Amersham Biosciences). The 
membranes were incubated with primary antibodies 
overnight at 4°C (DNMT1 [Abeam], MAT2A [Abeam], 
GAPDH [Santa Cruz]). Antibody binding was revealed 
by incubation with horseradish peroxidase-conjugated 
secondary antibodies followed by ECL Plus immuno- 
blotting detection system (Amersham Biosciences). 
Chemiluminescence was detected using a FluorChem 
HD2 camera with FluorChem software (Cell 
Biosciences). The signals were quantified using NIH 
Image J64 software and normalized relative to GAPDH. 

Statistical analysis. All experiments, except the 
Illumina® BeadChip determination of the DNA 
methylation pattern, included three biological 
replications, and statistically significant differences 
were determined by ANOVA (p < 0.05). For the 
Illumina HumanMethylation27 BeadChip, three 
replicates per population doubling level were studied 
and significant difference in site-specific methylation 
levels were determined as described above. 



1 www.funnet.info 

2 www.genecards.org 
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